In [1]:
import pandas as pd
import geopandas as gpd
from shapely.geometry import Point
import re
from tqdm.auto import tqdm
import contextily as cx
cx.set_cache_dir("tile_cache")
df = pd.read_excel("Pollution analyses Mapua.xlsx")
df.rename(columns={"Aluminuim": "Aluminium"})
df
Out[1]:
Sediment Location Total DDT Isomers Dieldrin Aluminuim Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
0 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.0075 0.001 16518.912655 1052.643599 42.463242 49.985737 273.474991 8.104016 52.778163 8.868234 35.506181 5.43789 0.991537 0.089184 0.016251 0.154799 8.519661 0.577402
1 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.0085 0.001 13182.681372 842.416269 36.994473 45.335967 191.183969 5.621974 35.468645 6.189135 25.626771 4.542981 0.980962 0.262338 0.01262 0.11341 5.90972 0.568269
2 Site 3 41˚ 15.493'S, 173˚ 5.876'E NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
3 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.014 0.001 14347.472272 865.1912 38.924618 42.498292 245.347282 6.065827 36.549068 6.342976 29.559022 6.078583 1.004799 0.026049 0.012362 0.027869 6.134818 0.44031
4 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.006 0.001 11359.826166 398.808354 38.351538 35.11283 211.319577 5.625269 31.457117 5.143459 23.527645 12.577909 0.867042 0.039512 0.010175 0.016312 7.136776 0.511512
5 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.006 0.001 12756.337883 619.394774 50.671753 40.737497 375.535486 7.412915 41.57042 6.000838 27.053254 27.206014 1.024441 0.03496 0.012906 0.017003 11.635085 0.569886
6 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.006 0.001 15457.058054 899.500095 40.896228 38.923171 279.821146 7.684669 44.240518 8.256869 32.816672 11.471663 0.97767 0.077509 0.014367 0.018943 7.481523 0.573935
7 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.006 0.001 17890.662309 998.201546 45.404394 43.747367 427.318202 8.376312 46.478442 9.123958 35.343672 14.571893 1.142411 0.042442 0.017004 0.020126 8.038082 0.68719
8 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0075 0.001 17473.290669 905.622552 45.969994 49.265992 272.912776 7.716316 48.967096 8.753198 38.218362 7.303673 1.159622 0.069248 0.014839 0.021722 8.128869 0.700793
9 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.006 0.001 25411.662513 1709.858197 57.596293 63.408278 386.369163 10.613864 69.300728 13.314035 50.822242 9.160778 1.438946 0.037045 0.020781 0.028689 11.237761 0.705635
10 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
11 Titiko - mud snails NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
12 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.1119 0.0048 9278.516506 10438.545347 26.122999 27.650273 191.842282 6.302344 37.37869 75.247891 49.645651 13.887514 3.036958 0.176368 0.052238 0.289213 5.886132 0.307983
13 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.1739 0.0083 3045.746777 9779.968577 9.266928 10.447118 87.293664 2.184617 15.890481 55.944349 32.241964 13.957774 2.454535 0.175301 0.040901 0.37774 1.59775 0.163822
14 Site 3 41˚ 15.493'S, 173˚ 5.876'E 0.2742 0.031 1688.387959 8956.741669 7.28127 13.902505 168.717108 1.792196 8.471845 47.170575 27.372176 14.702333 2.223494 0.08714 0.02785 0.181421 1.095187 0.125211
15 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.4465 0.023 5834.559801 6688.965169 16.851306 17.802685 176.232594 4.467314 30.757799 37.809536 32.459531 13.655855 1.834363 0.20797 0.076393 0.216377 3.078978 0.251346
16 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.0452 0.0026 776.544571 4906.246349 2.118771 2.044623 94.758462 0.811912 3.536308 19.422672 10.663063 6.717349 1.087014 0.148339 0.017945 0.107886 0.532221 0.061831
17 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0309 0.0015 1498.423285 7484.46063 4.984812 7.125436 114.022396 1.854082 7.87016 50.283271 25.967419 14.104914 2.469336 0.28798 0.040168 0.312954 1.141666 0.112747
18 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.0251 0.0015 2461.750292 13737.544243 6.892372 7.129266 57.248232 1.750822 9.710237 62.373216 37.50133 25.564629 2.140799 0.451634 0.192853 0.09457 1.821799 0.203198
19 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.0276 0.0005 66.04758 12226.537642 0.39743 0.301674 15.563173 0.744519 5.384783 64.524262 29.971684 10.063169 2.498558 0.550304 0.033397 0.163524 0.352246 0.046921
20 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0774 0.0015 2677.861843 7411.325007 6.881487 7.085678 63.300903 1.397971 8.251588 37.786295 23.023848 10.121427 1.267161 0.047833 0.026337 0.141486 1.747061 0.108743
21 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.0431 0.0005 4030.26235 11537.183933 9.555477 9.870691 84.961609 1.806037 11.641405 48.033265 26.564529 9.444581 2.063936 0.179229 0.025591 0.21965 1.855676 0.136361
22 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
23 Tuangi - cockles NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
24 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.0005 0.0005 787.639259 18087.464299 3.173642 3.555638 77.752051 2.079862 19.787828 6.039858 60.809418 30.604716 4.32808 0.091581 0.163688 0.088335 0.413524 0.2013
25 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.0005 0.0005 710.746254 17181.648407 2.634034 4.389376 20.956301 1.691787 15.865414 4.544055 49.425053 41.10754 3.770683 0.078458 0.155236 0.102417 0.437715 0.297476
26 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0005 0.0005 1166.470055 19532.038518 3.126026 3.65192 31.545677 1.952262 20.830599 5.549428 64.669416 23.676406 3.882477 0.184304 0.146906 0.077301 0.588409 0.151768
27 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
28 Patiki - Flounder Liver NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
29 Site 11 41˚ 15.680'S, 173˚ 5.875'E 0.0317 0.0014 42.706757 12777.79052 0.46972 0.184253 10.980888 1.146335 0.522057 26.782493 136.388133 24.703859 5.558432 0.183492 0.377747 0.376774 0.073363 0.012373
30 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
31 Patiki - Flounder Flesh NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
32 Site 11 41˚ 15.680'S, 173˚ 5.875'E 0.0042 0.0005 47.969064 10468.807791 0.173358 0.225523 1.891603 0.045363 0.225564 0.575889 37.250104 8.888757 1.288579 0.036508 0.00185 0.235958 0.035522 0.001545
33 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
34 Patiki - Flounder Flesh NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
35 NaN NaN NaN NaN 832.049226 9467.532397 2.155475 2.051407 48.205248 0.416026 2.524552 1.61799 81.990261 2.942643 1.165397 0.044176 0.007327 0.141226 0.308808 0.033755
36 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
37 Kutai - green-lipped mussel NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt Nickel Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium
38 Site 12 41˚ 15.704'S, 173˚ 5.959'E 0.003 0.0005 2382.523423 15487.941993 6.305763 6.132179 44.890628 1.771959 8.265882 5.052864 63.48187 13.710588 3.098001 0.040551 0.13091 0.157898 1.244749 0.14311
In [2]:
def parse_latlong(latlong):
    if type(latlong) is str:
        lat, long = re.findall(r"(\d+)˚ ([\d.]+)'(.)", latlong)
        lat_deg, lat_min, lat_dir = lat
        lon_deg, lon_min, lon_dir = long
        latitude = (float(lat_deg) + float(lat_min) / 60) * (-1 if lat_dir == 'S' else 1)
        longitude = (float(lon_deg) + float(lon_min) / 60) * (-1 if lon_dir == 'W' else 1)
        return Point(longitude, latitude)
str(parse_latlong("41˚ 15.624'S, 173˚ 5.892'E"))
Out[2]:
'POINT (173.0982 -41.2604)'
In [3]:
df["geometry"] = df.Location.apply(parse_latlong)
df = gpd.GeoDataFrame(df)
df.crs = 4326
df
Out[3]:
Sediment Location Total DDT Isomers Dieldrin Aluminuim Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium geometry
0 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.0075 0.001 16518.912655 1052.643599 42.463242 49.985737 273.474991 8.104016 ... 8.868234 35.506181 5.43789 0.991537 0.089184 0.016251 0.154799 8.519661 0.577402 POINT (173.09820 -41.26040)
1 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.0085 0.001 13182.681372 842.416269 36.994473 45.335967 191.183969 5.621974 ... 6.189135 25.626771 4.542981 0.980962 0.262338 0.01262 0.11341 5.90972 0.568269 POINT (173.09822 -41.25953)
2 Site 3 41˚ 15.493'S, 173˚ 5.876'E NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN POINT (173.09793 -41.25822)
3 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.014 0.001 14347.472272 865.1912 38.924618 42.498292 245.347282 6.065827 ... 6.342976 29.559022 6.078583 1.004799 0.026049 0.012362 0.027869 6.134818 0.44031 POINT (173.09652 -41.25838)
4 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.006 0.001 11359.826166 398.808354 38.351538 35.11283 211.319577 5.625269 ... 5.143459 23.527645 12.577909 0.867042 0.039512 0.010175 0.016312 7.136776 0.511512 POINT (173.09455 -41.25857)
5 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.006 0.001 12756.337883 619.394774 50.671753 40.737497 375.535486 7.412915 ... 6.000838 27.053254 27.206014 1.024441 0.03496 0.012906 0.017003 11.635085 0.569886 POINT (173.09288 -41.25860)
6 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.006 0.001 15457.058054 899.500095 40.896228 38.923171 279.821146 7.684669 ... 8.256869 32.816672 11.471663 0.97767 0.077509 0.014367 0.018943 7.481523 0.573935 POINT (173.08848 -41.25883)
7 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.006 0.001 17890.662309 998.201546 45.404394 43.747367 427.318202 8.376312 ... 9.123958 35.343672 14.571893 1.142411 0.042442 0.017004 0.020126 8.038082 0.68719 POINT (173.08700 -41.25658)
8 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0075 0.001 17473.290669 905.622552 45.969994 49.265992 272.912776 7.716316 ... 8.753198 38.218362 7.303673 1.159622 0.069248 0.014839 0.021722 8.128869 0.700793 POINT (173.09685 -41.25855)
9 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.006 0.001 25411.662513 1709.858197 57.596293 63.408278 386.369163 10.613864 ... 13.314035 50.822242 9.160778 1.438946 0.037045 0.020781 0.028689 11.237761 0.705635 POINT (173.09020 -41.25885)
10 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
11 Titiko - mud snails NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
12 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.1119 0.0048 9278.516506 10438.545347 26.122999 27.650273 191.842282 6.302344 ... 75.247891 49.645651 13.887514 3.036958 0.176368 0.052238 0.289213 5.886132 0.307983 POINT (173.09820 -41.26040)
13 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.1739 0.0083 3045.746777 9779.968577 9.266928 10.447118 87.293664 2.184617 ... 55.944349 32.241964 13.957774 2.454535 0.175301 0.040901 0.37774 1.59775 0.163822 POINT (173.09822 -41.25953)
14 Site 3 41˚ 15.493'S, 173˚ 5.876'E 0.2742 0.031 1688.387959 8956.741669 7.28127 13.902505 168.717108 1.792196 ... 47.170575 27.372176 14.702333 2.223494 0.08714 0.02785 0.181421 1.095187 0.125211 POINT (173.09793 -41.25822)
15 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.4465 0.023 5834.559801 6688.965169 16.851306 17.802685 176.232594 4.467314 ... 37.809536 32.459531 13.655855 1.834363 0.20797 0.076393 0.216377 3.078978 0.251346 POINT (173.09652 -41.25838)
16 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.0452 0.0026 776.544571 4906.246349 2.118771 2.044623 94.758462 0.811912 ... 19.422672 10.663063 6.717349 1.087014 0.148339 0.017945 0.107886 0.532221 0.061831 POINT (173.09455 -41.25857)
17 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0309 0.0015 1498.423285 7484.46063 4.984812 7.125436 114.022396 1.854082 ... 50.283271 25.967419 14.104914 2.469336 0.28798 0.040168 0.312954 1.141666 0.112747 POINT (173.09288 -41.25860)
18 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.0251 0.0015 2461.750292 13737.544243 6.892372 7.129266 57.248232 1.750822 ... 62.373216 37.50133 25.564629 2.140799 0.451634 0.192853 0.09457 1.821799 0.203198 POINT (173.08848 -41.25883)
19 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.0276 0.0005 66.04758 12226.537642 0.39743 0.301674 15.563173 0.744519 ... 64.524262 29.971684 10.063169 2.498558 0.550304 0.033397 0.163524 0.352246 0.046921 POINT (173.08700 -41.25658)
20 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0774 0.0015 2677.861843 7411.325007 6.881487 7.085678 63.300903 1.397971 ... 37.786295 23.023848 10.121427 1.267161 0.047833 0.026337 0.141486 1.747061 0.108743 POINT (173.09685 -41.25855)
21 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.0431 0.0005 4030.26235 11537.183933 9.555477 9.870691 84.961609 1.806037 ... 48.033265 26.564529 9.444581 2.063936 0.179229 0.025591 0.21965 1.855676 0.136361 POINT (173.09020 -41.25885)
22 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
23 Tuangi - cockles NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
24 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.0005 0.0005 787.639259 18087.464299 3.173642 3.555638 77.752051 2.079862 ... 6.039858 60.809418 30.604716 4.32808 0.091581 0.163688 0.088335 0.413524 0.2013 POINT (173.09820 -41.26040)
25 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.0005 0.0005 710.746254 17181.648407 2.634034 4.389376 20.956301 1.691787 ... 4.544055 49.425053 41.10754 3.770683 0.078458 0.155236 0.102417 0.437715 0.297476 POINT (173.09822 -41.25953)
26 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0005 0.0005 1166.470055 19532.038518 3.126026 3.65192 31.545677 1.952262 ... 5.549428 64.669416 23.676406 3.882477 0.184304 0.146906 0.077301 0.588409 0.151768 POINT (173.09288 -41.25860)
27 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
28 Patiki - Flounder Liver NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
29 Site 11 41˚ 15.680'S, 173˚ 5.875'E 0.0317 0.0014 42.706757 12777.79052 0.46972 0.184253 10.980888 1.146335 ... 26.782493 136.388133 24.703859 5.558432 0.183492 0.377747 0.376774 0.073363 0.012373 POINT (173.09792 -41.26133)
30 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
31 Patiki - Flounder Flesh NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
32 Site 11 41˚ 15.680'S, 173˚ 5.875'E 0.0042 0.0005 47.969064 10468.807791 0.173358 0.225523 1.891603 0.045363 ... 0.575889 37.250104 8.888757 1.288579 0.036508 0.00185 0.235958 0.035522 0.001545 POINT (173.09792 -41.26133)
33 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
34 Patiki - Flounder Flesh NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
35 NaN NaN NaN NaN 832.049226 9467.532397 2.155475 2.051407 48.205248 0.416026 ... 1.61799 81.990261 2.942643 1.165397 0.044176 0.007327 0.141226 0.308808 0.033755 None
36 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN None
37 Kutai - green-lipped mussel NaN Total DDT Isomers Dieldrin Aluminium Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium None
38 Site 12 41˚ 15.704'S, 173˚ 5.959'E 0.003 0.0005 2382.523423 15487.941993 6.305763 6.132179 44.890628 1.771959 ... 5.052864 63.48187 13.710588 3.098001 0.040551 0.13091 0.157898 1.244749 0.14311 POINT (173.09932 -41.26173)

39 rows × 21 columns

In [4]:
float_cols = df.columns[2:-1]
float_cols
Out[4]:
Index(['Total DDT Isomers', 'Dieldrin', 'Aluminuim', 'Sulphur', 'Vanadium',
       'Chromium', 'Manganese', 'Cobalt', 'Nickel', 'Copper', 'Zinc',
       'Arsenic', 'Selenium', 'Silver', 'Cadmium', 'Mercury', 'Lead',
       'Uranium'],
      dtype='object')
In [5]:
sediment = df.iloc[:10].copy()
sediment[float_cols] = sediment[float_cols].astype(float)
sediment
Out[5]:
Sediment Location Total DDT Isomers Dieldrin Aluminuim Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium geometry
0 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.0075 0.001 16518.912655 1052.643599 42.463242 49.985737 273.474991 8.104016 ... 8.868234 35.506181 5.437890 0.991537 0.089184 0.016251 0.154799 8.519661 0.577402 POINT (173.09820 -41.26040)
1 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.0085 0.001 13182.681372 842.416269 36.994473 45.335967 191.183969 5.621974 ... 6.189135 25.626771 4.542981 0.980962 0.262338 0.012620 0.113410 5.909720 0.568269 POINT (173.09822 -41.25953)
2 Site 3 41˚ 15.493'S, 173˚ 5.876'E NaN NaN NaN NaN NaN NaN NaN NaN ... NaN NaN NaN NaN NaN NaN NaN NaN NaN POINT (173.09793 -41.25822)
3 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.0140 0.001 14347.472272 865.191200 38.924618 42.498292 245.347282 6.065827 ... 6.342976 29.559022 6.078583 1.004799 0.026049 0.012362 0.027869 6.134818 0.440310 POINT (173.09652 -41.25838)
4 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.0060 0.001 11359.826166 398.808354 38.351538 35.112830 211.319577 5.625269 ... 5.143459 23.527645 12.577909 0.867042 0.039512 0.010175 0.016312 7.136776 0.511512 POINT (173.09455 -41.25857)
5 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0060 0.001 12756.337883 619.394774 50.671753 40.737497 375.535486 7.412915 ... 6.000838 27.053254 27.206014 1.024441 0.034960 0.012906 0.017003 11.635085 0.569886 POINT (173.09288 -41.25860)
6 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.0060 0.001 15457.058054 899.500095 40.896228 38.923171 279.821146 7.684669 ... 8.256869 32.816672 11.471663 0.977670 0.077509 0.014367 0.018943 7.481523 0.573935 POINT (173.08848 -41.25883)
7 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.0060 0.001 17890.662309 998.201546 45.404394 43.747367 427.318202 8.376312 ... 9.123958 35.343672 14.571893 1.142411 0.042442 0.017004 0.020126 8.038082 0.687190 POINT (173.08700 -41.25658)
8 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0075 0.001 17473.290669 905.622552 45.969994 49.265992 272.912776 7.716316 ... 8.753198 38.218362 7.303673 1.159622 0.069248 0.014839 0.021722 8.128869 0.700793 POINT (173.09685 -41.25855)
9 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.0060 0.001 25411.662513 1709.858197 57.596293 63.408278 386.369163 10.613864 ... 13.314035 50.822242 9.160778 1.438946 0.037045 0.020781 0.028689 11.237761 0.705635 POINT (173.09020 -41.25885)

10 rows × 21 columns

In [6]:
for col in tqdm(float_cols):
  ax = sediment.plot(col, figsize=(10,10), legend=True, legend_kwds={"orientation": "horizontal", "pad": 0.05}, cmap="Reds", markersize = 100, edgecolor = "black")
  ax.set_title(f"Sediment {col} in Mapua")
  ax.ticklabel_format(useOffset=False, style='plain')
  xmin, ymin, xmax, ymax = sediment.total_bounds
  pad = .002
  ax.set_xlim(xmin-pad, xmax+pad)
  ax.set_ylim(ymin-pad, ymax+pad)
  cx.add_basemap(ax, source=cx.providers.Esri.WorldImagery, crs=sediment.crs, zoom=18)
  file_name = f"Sediment {col} in Mapua.png"
  ax.figure.savefig(file_name, bbox_inches="tight")

  0%|          | 0/18 [00:00<?, ?it/s]
In [7]:
titiko = df.iloc[12:22].copy()
titiko[float_cols] = titiko[float_cols].astype(float)
titiko
Out[7]:
Sediment Location Total DDT Isomers Dieldrin Aluminuim Sulphur Vanadium Chromium Manganese Cobalt ... Copper Zinc Arsenic Selenium Silver Cadmium Mercury Lead Uranium geometry
12 Site 1 41˚ 15.624'S, 173˚ 5.892'E 0.1119 0.0048 9278.516506 10438.545347 26.122999 27.650273 191.842282 6.302344 ... 75.247891 49.645651 13.887514 3.036958 0.176368 0.052238 0.289213 5.886132 0.307983 POINT (173.09820 -41.26040)
13 Site 2 41˚ 15.572'S, 173˚ 5.893'E 0.1739 0.0083 3045.746777 9779.968577 9.266928 10.447118 87.293664 2.184617 ... 55.944349 32.241964 13.957774 2.454535 0.175301 0.040901 0.377740 1.597750 0.163822 POINT (173.09822 -41.25953)
14 Site 3 41˚ 15.493'S, 173˚ 5.876'E 0.2742 0.0310 1688.387959 8956.741669 7.281270 13.902505 168.717108 1.792196 ... 47.170575 27.372176 14.702333 2.223494 0.087140 0.027850 0.181421 1.095187 0.125211 POINT (173.09793 -41.25822)
15 Site 4 41˚ 15.503'S, 173˚ 5.791'E 0.4465 0.0230 5834.559801 6688.965169 16.851306 17.802685 176.232594 4.467314 ... 37.809536 32.459531 13.655855 1.834363 0.207970 0.076393 0.216377 3.078978 0.251346 POINT (173.09652 -41.25838)
16 Site 5 41˚ 15.514'S, 173˚ 5.673'E 0.0452 0.0026 776.544571 4906.246349 2.118771 2.044623 94.758462 0.811912 ... 19.422672 10.663063 6.717349 1.087014 0.148339 0.017945 0.107886 0.532221 0.061831 POINT (173.09455 -41.25857)
17 Site 6 41˚ 15.516'S, 173˚ 5.573'E 0.0309 0.0015 1498.423285 7484.460630 4.984812 7.125436 114.022396 1.854082 ... 50.283271 25.967419 14.104914 2.469336 0.287980 0.040168 0.312954 1.141666 0.112747 POINT (173.09288 -41.25860)
18 Site 7 41˚ 15.530'S, 173˚ 5.309'E 0.0251 0.0015 2461.750292 13737.544243 6.892372 7.129266 57.248232 1.750822 ... 62.373216 37.501330 25.564629 2.140799 0.451634 0.192853 0.094570 1.821799 0.203198 POINT (173.08848 -41.25883)
19 Site 8 41˚ 15.395'S, 173˚ 5.220'E 0.0276 0.0005 66.047580 12226.537642 0.397430 0.301674 15.563173 0.744519 ... 64.524262 29.971684 10.063169 2.498558 0.550304 0.033397 0.163524 0.352246 0.046921 POINT (173.08700 -41.25658)
20 Site 9 41˚ 15.513'S, 173˚ 5.811'E 0.0774 0.0015 2677.861843 7411.325007 6.881487 7.085678 63.300903 1.397971 ... 37.786295 23.023848 10.121427 1.267161 0.047833 0.026337 0.141486 1.747061 0.108743 POINT (173.09685 -41.25855)
21 Site 10 41˚ 15.531'S, 173˚ 5.412'E 0.0431 0.0005 4030.262350 11537.183933 9.555477 9.870691 84.961609 1.806037 ... 48.033265 26.564529 9.444581 2.063936 0.179229 0.025591 0.219650 1.855676 0.136361 POINT (173.09020 -41.25885)

10 rows × 21 columns

In [8]:
for col in tqdm(float_cols):
  ax = titiko.plot(col, figsize=(10,10), legend=True, legend_kwds={"orientation": "horizontal", "pad": 0.05}, cmap="Reds", markersize = 100, edgecolor = "black")
  ax.set_title(f"Tītiko (mud-flat snails) {col} in Mapua")
  ax.ticklabel_format(useOffset=False, style='plain')
  xmin, ymin, xmax, ymax = titiko.total_bounds
  pad = .002
  ax.set_xlim(xmin-pad, xmax+pad)
  ax.set_ylim(ymin-pad, ymax+pad)
  cx.add_basemap(ax, source=cx.providers.Esri.WorldImagery, crs=sediment.crs, zoom=18)
  file_name = f"Tītiko (mud-flat snails) {col} in Mapua.png"
  ax.figure.savefig(file_name, bbox_inches="tight")

  0%|          | 0/18 [00:00<?, ?it/s]